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Spherically symmetric inhomogeneous dust collapse has been studied in higher dimensional space- 
time and the factors responsible for the appearance of a naked singularity are analyzed in the region 
close to the centre for the marginally bound case. It is clearly demonstrated that in the former case 
naked singularities do not appear in the space-time having more than five dimension, which appears 
to a strong result. The non-marginally bound collapse is also examined in five dimensions and the 
role of shear in developing naked singularities in this space-time is discussed in details. The five 
dimensional space-time is chosen in the later case because we have exact solution in closed form 
only in five dimension and not in any other case. 
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I. INTRODUCTION 

The Cosmic Censorship Conjecture [1,2] is yet one of the unresolved problems in General 
Relativity. It states that the space-time singularity produced by gravitational collapse must be 
covered by the horizon. However the singularity theorems as such do not state anything about the 
visibility of the singularity to an outside observer. In fact the conjecture has not yet any precise 
mathematical proof. Several models related to the gravitational collapse of matter has so far been 
constructed where one encounters a naked singularity [3-8]. 

It has recently been pointed out by Joshi et al [9] that the physical feature which is re¬ 
sponsible for the formation of naked singularity is nothing but the presence of shear. It is 
the shear developing in the gravitational collapse, which delays the formation of the apparent 
horizon so that the communication is possible from the very strong gravity region to observers 
situated outside. Joshi et al have analyzed in details how the presence of shear determines the 
growth and evolution of an inhomogeneous dust distribution represented by four dimensional 
Tolman-Bondi metric and have attempted to clarify the nature of singularities as the final outcome. 

The objective of this paper is to fully investigate the situation in the background of higher 
dimensional space-time with arbitrary dimensions for marginally bound collapse and in five 
dimensional space-time for non-marginally bound case. The reason for confining our discussions to 
five dimensions only in the second case is that the closed form solutions for non-marginally bound 
collapse are available only in 5 D space-time. These discussions are of adequate relevance in the 
context of recent interest generated in the study of gravitational collapse in higher dimensional 
space-time [10-16]. 

We have shown in the present paper that for marginally bound collapse that is for / = 0 the 
naked singularity may appear only when the space-time has dimensions upto five which appears 
to a strong result. In more than five dimensions t a h < to that is the apparent horizon at any 
point in the region forms earlier than the central shell focusing singularity, which indicates that 
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the shell focusing singularity first appearing at r = 0 remains hidden behind the apparent horizon 
and so gives rise to a black hole always. 

The non-marginally bound collapse of a 5 D inhomogeneous dust reveals that there may occur 
both the black hole and the naked singularities depending on the initial data. The five dimensional 
space-time is chosen in the later case because we have exact solution in closed form only in five 
dimension and not in any other case. 


II. HIGHER DIMENSIONAL INHOMOGENEOUS DUST 

The higher dimensional Tolman-Bondi type metric is given by 


ds 2 = e v dt 2 — e x dr 2 — R 2 dVl 2 n (1) 

where u, A, R are functions of the radial co-ordinate r and time t and dtt^ represents the metric 
on the n-sphere. Since we assume the matter in the form of dust, the motion of particles will be 
geodesic allowing us to write 


= 1 


( 2 ) 


Using comoving co-ordinates one can in view of the field equations [11] arrive at the following 
relations in (n + 2) dimensional space-time 


and 


R' 2 

1 + f{r) 


( 3 ) 


R 2 


f(r) 


F(r) 
f ?"- 1 ’ 


( 4 ) 


where the function f(r) classifies the space-time [2] as bound, marginally bound and unbound 
depending on the range of its values which are respectively 


f(r) < 0, /(r) = 0, f(r) > 0. (5) 

The function F(r) can be interpreted as the mass function which is related with the mass 
contained within the comoving radius r. 

The above models are characterized by the initial data specified on the initial hypersurface t = ti, 
from which the collapse develops. As it is possible to make an arbitrary relabeling of spherical 
dust shells by r —> g(r), without loss of generality we fix the labeling by the choice R(U, r ) = r, so 
that initial density distribution is given by 


Pi(r) = 


nF' 
2r” 


( 6 ) 


Now the integration of the equation (4) gives us the exact solution for I? as a function of r and 
t. For the case / = 0 the integration is straightforward to yield the solution 


R = 


+ 1 ^ n+1 


F n + 1 (■t 


S 


ty 


( 7 ) 


where t s is a function of the radial co-ordinate r and t = t s is the instant of shell focusing 
singularity occurring at r, that is R(t s (r),r) = 0. It indicates that the shell focusing singularity 
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occurs at different r at different epochs. So the collapse in this case is not simultaneous in comoving 
co-ordinates and is in variance with that in the homogeneous dust model [6,7,12]. The second case 
is for / ^ 0. One can not, however, obtain the explicit solution of (4) in closed form except for 
n = 3 that is in five dimensional space-time. Even in the simplest case of 4 dimensions the solution 
for non-zero f(r) is obtained in the parametric form only. 


III. MARGINALLY BOUND CASE (/ = 0) IN (n + 2 )D SPACE-TIME 

We observe that the normal vector to the boundary of any hypersphere R — Rq = 0 is given by 

In = (R, R', 0,0, 0,.). In order that is a null vector and the boundary of this hypersphere is 

a null surface we must have 


R(t ah ,r) = [F(r)}^~ (8) 

where t = t a h{r) is the equation of apparent horizon which marks the boundary of the trapped 
region. If the apparent horizon develops earlier than the time of singularity formation, the event 
horizon can fully cover the singularity, which then may be said to be hidden within a black hole. 
Combining (7) and (8) along with the initial condition R(ti,r) = r we arrive at the relation 


(tah 


ti) 


2 

n + 1 


n+1 - 

r 2 F 


1/2 


2 

n + 1 


1 

p— 


(9) 


Since our study is restricted to the region near r = 0, the mass function F(r) should vanishes 
exactly at r = 0 and it is possible to express F(r) as a polynomial function of r near the origin. 
One should at the same time keep in mind that initially the central density is regular. Further from 
physical considerations one may argue that p' (r) should vanish exactly at r = 0, but is negative 
in the neighbouring region. All these consideration lead to the following expressions for F(r) and 
the initial energy density p/, 


F(r) = [F 0 r n+1 + F 2 r n+3 + F 3 r”+ 4 

and 


( 10 ) 


Tl r> 

Pi( r ) = 2 [( n + X ) F o + ( n + 3)^2 r +.] (11) 

Obviously the initial central density is given by p c = ”( rt + 1 ) p 0 . Since p'(r) < 0 in the region 
r « 0 one must have F 2 < 0. 

Now in view of (7) one can write 


r (n+l)/2 


I'S ti 


(n + 1) [F 0 r n + 1 + F 2 r n + 3 + ...] 1/2 


( 12 ) 


If we now denote t s = to as the instant of first shell focusing singularity to occur at the centre 
r = 0, the relation (12) can be used to obtain the following relation 


to — U — 


(n + W 

which in combination with the relation (9) yields 


1/2 ’ 


(13) 
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The relation (14) is valid in general (n + 2) dimensional space-time and reduces to the 4 dimen¬ 
sional expression (n = 2) in the region near the centre as given in the paper of Joshi et al [9]. In 
5 D space-time the naked singularity appears only if F 2 ^ 0 and further F 2 | > 2 Fq . Otherwise for 
either F 2 = 0 or |F 2 | < 2 Fq there must occur a black hole, because in that case t a h < to that is, 
the horizon appears earlier than the shell focusing singularity at r = 0. Only in a very special case 
of |F 2 | = 2 Fq the occurrence of naked singularity or black hole will depend on the co-efficients of 
higher powers of r. For general (n + 2) dimensions the close examination of (14) reveals that if 
there is to exist a naked singularity we must have 


n + 1 
n — 1 


> 2 , 


(15) 


which means n < 3. It is interesting to note that has the value 2 when n = 3 and 

then it decreases monotonically with the increasing number of dimensions. So for space-time 

i n _|_ i 

with dimension larger than five the term r"- 1 dominates leading to a black hole. Our 

conclusion is that for the space-time having larger than five dimensions the existence of the 
naked singularity is prohibited or in other words the shell focusing singularity is fully covered 
by the horizon in such cases. However, if we relax the restriction namely, p'(r) = 0 at r = 0 
then F\ will be non-zero and this will leads to the possibility of naked singularity in any dimension. 


Calculation of shear near r = 0: 


We still consider the marginally bound case / = 0. The shear scalar in the (n + 2) dimensional 
spherically symmetric dust metric so far discussed in this section may be estimated by the factor 
[9] a = y/ 2(n w | _ 1) • Using (3) and (4) one gets in turn 

&_R\ = / n [RF 1 - (n + 1)#F] 

R' R) V 8 (n + 1) F 1 /2 j R(n+i)/2 i? / 1 > 

From (7) 


2 (n + 1) 


R' = 




t) "+ 1 + 


n +1 


F"+ 1 (t s 


t)- 


-l 


t 


/ 
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where F,F',t s etc.’s are all functions of the comoving radial co-ordinate r. So as t approaches 
t s , R approaches zero and effectively cr 2 approaches (n+1) 2 R ^ +1 , which in turn becomes infinitely 
large. It is therefore evident that at each r, the magnitude of shear explodes as the energy density 
explodes (R = 0) and this occurs at different instants at different spherical shells of different radial 
co-ordinates including the centre (r = 0). 


It is now possible in view of the choice R(U, r) = r to express the initial shear cr,; to the form 


<Ji = 


'[rF' - (n + 1)F] 


8(n + 1) F 1 / 2 r("+ 1 )/ 2 

But using (7) along with the initial condition R(ti, r ) = r, we have a relation like 


(17) 


n+1 r. 

r~ -R- 


t-U = F 


1/2 


(18) 
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which in turn being used in (16) yields the expression for shear scalar a(t, r ) very near the centre 
(r « 0) in the following form after some manipulations: 



so the initial shear at t = U is given by 



We find that at the centre (r = 0) the initial shear vanishes. It is interesting to note that the 
dependence on r of the initial shear does not depend on the number of dimensions of the space- 
time. In fact the expression (20) coincides exactly with the value of the initial shear ex,; calculated 
earlier by Joshi et al [9] for 4 D space-time (n = 2). In view of what has been discussed so far 
the expression (20) reveals that the existence of the naked singularity is directly related with the 
non-vanishing shear in four and five dimensions. But in larger dimensions the central singularity 
seems to be covered by the appearance of apparent horizons. One should mention here that the 
statement of the paper of Joshi et al [9] that the shear decreases in course of time and finally 
vanishes at t = to is not true. In fact the shear in view of (19) increases and goes to infinity at 
this epoch, which is expected. 

IV. NON-MARGINALLY BOUND FIVE DIMENSIONAL SPACE-TIME (/ ^ 0) 

We already know that for / ^ 0 it is possible to obtain the solution of (4) in closed form only 
in 5 D space-time that is for n = 3. In other cases even for the simplest case of a four dimensional 
manifold the solution is available only in parametric form. In 5 D space-time one of the exact 
solution of the equation (4) is given by 


R 2 = [f(t s - t) 2 + 2F 1 / 2 (t s - i)] (21) 

For / = 0 the solution (21) is coincident with the solution (7) when we put n = 3. By the same 
arguments put forward in the previous section it is possible to express F(r) as 

F(r) = F 0 r 4 + F 2 r 6 + F 3 r 7 + . (22) 

and the initial density distribution pi{r) as 

Pi{r) = p c + p 2 r 2 + p 3 r 3 + ., (23) 

where and p c = 6 Fq, p 2 = 9F 2 , so on. 


Here also following the earlier reasoning F 2 < 0. The only new input in the present case is the 
function f(r), which may be expressed as a power series in r near the centre r = 0. 

We assume 

f(r) = for + fir 2 + f 2 r 3 + . (24) 

The form (24) is chosen because f(r) vanishes as r —> 0, which is demanded by the regularity 
condition at r = 0 [17]. 
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Now when t = ti one can write from (21) 


(F + /r 2 ) 1//2 F 1 / 2 

7 T 


( 25 ) 


The finite non zero magnitude of the left hand side of (25) demands fo = 0 and fi / 0 as 
minimum requirements. The expansion of (22) near r = 0 yields 




(*0 + A) 1/2 - Fl' 2 , f 2 


h 


+ 


P 


2F^ 2 (Fq + fi) 1 / 2 — (2Fq + fi) 
2 (F 0 + A)i/2 


/I 


2pV7T+Fb 


+ 


(P-fif3)W.h + F 0 -VK) f 2 , 4(/ 1 + f 0 )(/ 3 + f 2 )-/. 


/i 3 2/iVTb 8/i(/i + To) 3 ^ 2 

At the centre r = 0 the relation (26) reduces exactly to 


(26) 


to ~ ti 


(F 0 + fi) 1 ' 2 - -Fq 1/2 

h 


(27) 


Now at t = t a h one must have R 2 


F as justified earlier. Hence by solving (21) we get 


t a h(r) - t s (r) 


pl/2 _ p 1 / 2 ^ + jy/2 

7 


(28) 


The form (28) is chosen because it satisfies the consistency relation in the limit / ss 0. Remem¬ 
bering that t s (r = 0) = to we obtain in the region very close to the centre 


tah(r) = t 0 + 


ll 

P 


2F 0 1/2 (F 0 + A) 1 / 2 - (2F 0 + fi) 

2 ( Fo + / 1 ) 1 /2 


P 


2 PVJT+Fo 


fiF 0 + F 2 

vPM 


, (/I - fihXVKTFb - v^) , 4(/r + F 0 )(f 3 + F 2 ) - /|] , 

+ p + s/i (./i • /-o) :i /2 r . ( - J) 

In the above relation if / 2 < 0 the second term on the R.H.S is positive, which can easily be 
verified. Hence it is possible to conclude that t a h > to near r = 0. This suggests that that the 
apparent horizon in the region close to the centre appears later than the shell focusing singularity 
at r = 0 and hence the naked singularity should exist in this case. On the other hand if / 2 > 0 
it should be a black hole. These conclusions are valid irrespective of f(r) be positive or negative 
except for the restriction Fo > |/i| in case /i < 0. However, no definite conclusion can be drawn 
in case f 2 = 0. 

Calculation of shear: 

The general expression for the shear in five dimensional space-time is given by 

Using the expansions of F(r) and /(r) from (22) and (24) along with the assumption fo = 0, 
which has already been justified earlier we arrive at the expression for the shear at any instant in 
the region close to r = 0. This expression is given after omitting a few intermediate steps by 


/' F' 

~R' + R 2 R' ~ R3 


4 F 2/ 

R 


(30) 
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a = v ! (Fo+/i) " 1/2 


/ 2 , 

— r + m 
JC 


Rrn | fm +1 

X3 A' 


+ 0(r m+1 ) 


(31) 


where to > 2 and A = X(i) = [1 + fi(t — t;) 2 — 2(F 0 + /i) 1 ^ 2 (i — i*)]. 


Since er depends on the time the shear develops in course of time at any r away from the centre. 
When / = 0 we get for the initial shear at t = U 


3 mF, 


32 2 




+ 0(r m+i ) 


with to > 2. This is the same expression as in 4 D case. 

It is evident that the initial shear vanishes when f 2 = 0, (F 2 + fy) = (F^ + /*) = . = 0 and 

hence even if the initial shear is zero the dust distribution may be inhomogeneous because the 
co-efficient F 2 may still be non-zero. This makes in view of (26) t s a function of the comoving 
radial co-ordinate r, so that the shell focusing singularity appears at different r at different 
instants. The nature of singularities appearing in marginally bound cases (/ = 0) is therefore 
clearly distinct from the present case of non-marginally bound case (/ 7 ^ 0). This fact, however, 
has not been taken care of by Joshi etal in their discussion in the 4 dimensional marginally bound 
space-time. 

Acknowledgement: 

The authors are grateful to the referee for his valuable comments and suggestions to improve 
the manuscript. They further thank the members of the Relativity and Cosmology Research 
Centre, Department of Physics, Jadavpur University for helpful discussion. One of the authors 
(U.D) is thankful to CSIR (Govt, of India) for awarding a Junior Research Fellowship. 

References: 

[1] R. Penrose, Riv. Nuovo Cimento 1 252 (1969); in General Relativity, an Einstein Cen¬ 
tenary Volume, edited by S.W. Hawking and W. Israel (Cambridge Univ. Press, Cambridge, 
1979). 

[2] P.S. Joshi, Global Aspects in Gravitation and Cosmology (Oxford Univ. Press, Oxford, 1993). 

[3] B. Waugh and K. Lake, Phys. Rev. D 38 1315 (1988). 

[4] J. P. S. Lemos, Phys. Lett. A 158 271 (1991); Phys. Rev. Lett. 68 1447 (1992). 

[5] P. S. Joshi and I. H. Dwivedi, Class. Quantum Grav. 16 41 (1999). 

[ 6 ] A. Ilha and J. P. S. Lemos, Phys. Rev. D 55 1788 (1997). 

[7] A. Ilha, A. Kleber and J. P. S. Lemos, J. Math. Phys. 40 3509 (1999). 

[ 8 ] J. F. V. Rocha and A. Wang, ibid 17 2589 (2000). 

[9] P.S. Joshi, N. Dadhich and R. Maartens, Phys. Rev. D 65 101501(i?)(2002). 

[10] D. M. Eardley and L. Smarr, Phys. Rev. D 19 2239 (1979); D. M. Eardley in Gravitation 
and Astrophysics, edited by B. Carter and J. B. Hartle (NATO Advanced Study Institute, Series 
B; Physics, Vol. 156) (Plenum Press, New York, 1986), pp. 229-235. 

[11] A. Benerjee, A. Sil and S. Chatterjee, Astrophys. J. 422 681 (1994); A. Sil and S. Chatterjee, 
Gen. Rel. Grav. 26 999 (1994). 

[12] S. Chatterjee, A. Banerjee and B. Bhui, Phys. Lett. A 149 91 (1990). 

[13] S. G. Ghosh and N. Dadhich, Phys. Rev. D 64 047501 (2001). 

[14] S. G. Ghosh and A. Beesham, Phys. Rev. D 64 124005 (2001). 

[15] S. G. Ghosh and A. Beesham, Class. Quantum Grav. 17 4959 (2000). 

[16] S. S. Deshingkar, P.S. Joshi and I. H. Dwivedi, Phys. Rev. D 59 044018 (1999). 

[17] T. Harada, H. Iguchi and K.I. Nakao, Phys. Rev. D 61 101502 (2000); Prog. Theor. Phys. 
107 449 (2002) . 




